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The set of Faddeev and Lippmann-Schwinger integral equations for three-body systems involving 
Coulomb interactions deduced from a "three-potential" picture are shown to be compact for all 
energies and a method of solution is given. 
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There are two well-known, practical but genuinely dif- 
ferent approaches to solve three-body scattering equa- 
tions of Faddeev type and which involve Coulomb-like 
interactions. In the first of these, the Alt-Grassberger- 
Sandhas (AGS) equations are solved by the "screening 
and renormalizing" technique. This approach and its 
applications are discussed in a recent review [jjj. The 
alternative approach is to solve the Faddeev-Noble dif- 
ferential equations. To do so, in configuration space, one 
needs asymptotic boundary conditions. This approach 
has also been reviewed recently [Q. Although these ap- 
proaches have been in use for some time, problems still 
exist with the renormalization procedure of method one 
and with the approximate boundary condition required 
in method two. Consequently, to date, only limited solu- 
tions are available below or above the breakup threshold. 

Herein I propose a new approach for solving the three- 
body Coulomb scattering problem. As usual, one pre- 
sumes that the quantum mechanical system evolves from 
a state governed by the asymptotic Hamiltonian to the 
physical state described by the total Hamiltonian. In the 
two-potential formalism of that, an intermediate Hamil- 
tonian must be defined and connection made first to the 
asymptotic Hamiltonian and then to the total Hamilto- 
nian. The two-potential picture is a direct consequence 
of Kato's chain rule for wave operators Q, which holds 
even for Coulomb interactions Q . In a "three-potential" 
picture || a three-body Coulomb scattering process can 
be viewed as three consecutive scattering processes by 
which the asymptotic channel Hamiltonian connects to 
the total one via two intermediate Hamiltonians. This 
process can be formulated in terms of a set of the Fad- 
deev and Lippmann- Schwinger integral equations. Nu- 
merically, these integral equations can be solved by us- 
ing a Coulomb-Sturmian space representation. Hereafter 
I designate both the representation and the method by 
the label CS. The (CS) method has been used before and 
it has worked very well for bound-state problems with re- 
pulsive H and with attractive Coulomb interactions. 
It has also been used to analyze p — d scattering at en- 
ergies below breakup || and the results agree very well 
with those obtained from other calculations §31. 



In this article, I seek to extend the scope of this novel 
(CS) method to energies above the breakup threshold. I 
shall show that the Faddeev and Lippmann-Schwinger 
integral equations deduced from a "three-potential" pic- 
ture by applying Kato's chain rule possess compact ker- 
nels and so have unique solutions to the three-body 
Coulomb problem for all energies. I will also show how 
one can calculate the Green's operator which contains all 
the asymptotically relevant terms for arbitrary complex 
energies. 

I assume that the subsystem interaction, v a , of two 
elements in a three-body system is a Coulomb-like one, 
which is then split into short-range and long-range terms 
as 



««(&») =«i' } (&) + *£ J (&) 



(I), 



(1) 



Here a denotes the subsystem and superscripts s and I in- 
dicate the short- and long-range attributes, respectively, 
with £ a and r\ a being the usual configuration space Jacobi 
coordinates. The splitting should be performed in such 
a way that Va does not support any bound states. The 
total Hamiltonian is denoted by H and the asymptotic 
channel Hamiltonian is defined by 



H, 



v a = H° 



(2) 



where h and H° are the two-body and the three-body 
kinetic energy operators respectively. The asymptotic 
state |$ Q ) is an eigenstate of H a , i.e. 



H a \$ a ) =E a \* a ), 



(3) 



where ({.aVal^a) = (ValXa) (^a\4>a) is a product of a con- 
tinuum state in coordinate r\ a and a bound or contin- 
uum state in the two-body subsystem £ Q . If <p a is a 
bound eigenstate, asymptotically there is a two-fragment 
channel. If it is a scattering state, the asymptote is a 
three-fragment channel. The scattering state, $ a , which 
evolves from the asymptotic state, is given by applying 
the Q^(H,H a ) M0ller operators onto the asymptotic 
state 



(4) 
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The M0ller operators are defined by the limit 

n {±) (H, H a )=s- limexp [\Ht] exp [-iH a (t)], (5) 

t — >^oo 

where H a (t) = H a t + A(t) with A(t) being the Dollard's 
distortion operator ||. 

In the spirit of Kato's chain rule, I introduce two in- 
termediate Hamiltonians, the channel long-range Hamil- 
tonian 



H$ =H a +vf +v«\ (6) 
and the channel distorted long-range Hamiltonian 



H a = H a 



,(') 



(7) 



The auxiliary potential — Ua\r) a ) is defined such 
that it does not support any bound states and has the 
asymptotic form ~ e a (ep + e 7 )/7y a as r\ a — > oo. In 
fact, Ua is an effective Coulomb interaction between the 
center of mass of the subsystem a (with charge ep + e 7 ) 
and the third particle (with charge e Q ). Now, accord- 
ing to Kato's chain rule the M0ller operator can be 
written in the form 

xfl^(H^,H a )Q^(H a ,H a ). (8) 

The last term, fl^'(H a ,H a ), which essentially de- 
scribes a two-body scattering in the Coulomb-like po- 
tential u a (r] a ), exists in Dollard's sense, and by applying 
it to the channel state gives 



\¥*>) = ni ± \H a ,H a )\* a ) = \x£'>)\<i> a ) 



(9) 



with (rj a \xa^) being a scattering solution in the 
Coulomb-like potential u„ . 

Now I will show that the middle term in (^) exists in 
the ordinary sense jlO| 

n {±) (H2 ] ,H a ) = s - limexp [iH a l) t] exp [-\H a t]. (10) 

t — >=foo 

The necessary condition for the existence of the strong 
limit in jic|) is that the potential 



U a = ff« - H a = v 



(0 




,(0 



(11) 



should decrease in configuration space more rapidly than 
does the Coulomb potential Q. Indeed, the auxiliary 
potential u„ has been constructed so that it cancels 
the long-range potentials vV + in U a for the two- 
fragment channel a in which the particle labeled with a 
is in continuum state, while the particles labeled by (3 
and 7 form a bound state. Also, in the three-fragment 
channel of H a , where the particles labeled by j3 and 7 are 
in continuum state, U a decreases with radius faster than 



does the Coulomb potential. The term U a would behave 
like a Coulomb potential in the two- fragment channels, 
(3 and 7, where particles a and 7 or a and (3 form bound 
clusters, but for there are no two-fragment channels 
(3 and 7. 

Applying n^(H^>,H a ) to the "free" state |$i ±} ) 
yields 



(12) 



Sandhas has shown in Ref. |12[, that if for a three-body 
system only one of the potentials supports bound states, 
no rearrangement channels are possible, and therefore 
a single Lippmann-Schwinger equation is completely 
sufficient to guarantee the uniqueness of the solution. 
The Hamiltonian Ha possesses this property, thus the 
Lippmann-Schwinger integral equation which is derived 
from (0), 

|$«(±)) = |$(±)) + c a (E a ± iO)XJ a \$>WW), (13) 

with G a (z) = {z — Hq,) -1 , provides unique solution for 
this auxiliary three-body system. Additionally, Eq. ( p^ ) 
possesses a compact kernel because in the asymptotically 
accessible region, channel a, G a is linked solely to the 
asymptotic behavior of |$a' )( ' ± ' ) }, and the source term, 
U a , is of shorter range than the Coulomb interaction. 
Thus, in this particular Lippmann-Schwinger equation 
U a can be approximated by finite rank terms. Note, 



that a similar equation holds for Ga = (z — H a ' 
G® =G a + G a U a G ( a l) . 



r(0 



(14) 



Finally, Q^\H,HcP), which also exists in the usual 
sense, leads to the Faddeev-Noble integral equations [^3) , 

\4 ±] ) = <WI^ K±) > + Gf{E a ± i0)^ } E |^)), 

(15) 

where a, (3, 7 form a cyclic permutation. Merkuriev has 
shown that the kernels of these equations are compact 
for all energies ]l4| ], and thus the potential operators 
can be approximated well by finite rank terms. So, the 
set of Lippmann-Schwinger and Faddeev-Noble integral 
equations derived from a "three-potential" picture for 
the three-body Coulomb scattering problem, are uniquely 
solvable and possess compact kernels for all energies. 

In this derivation a crucial point is that one Lippmann 
Schwinger equation of the type given in Eq. (13) suffices 
for a unique solution to be found. This is the case if 
one imposes a condition on that it should not sup- 
port bound states. This condition is satisfied if all the 
Coulomb interactions in the three-body system are repul- 
sive. With a three-nucleon system then, one would take 
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Va ^ as the nuclear potential and Va as the Coulomb in- 
teraction. The situation is more complicated if some of 
the Coulomb interactions are attractive. There are an 
infinite number of bound states associated with an at- 
tractive Coulomb potential; a fact closely related to its 
long-range character. Nevertheless, the procedure is ap- 
plicable to such systems whenever the energy is below 
the three-body breakup threshold. For such cases one 
Lippmann-Schwinger equation is still sufficient to find 
solutions. 

The "three-potential" integral equations have been 
solved previously by using a CS representation 
With n and I being the radial and orbital angular mo- 
mentum quantum numbers respectively, the CS func- 
tions \nl) form a biorthonormal discrete basis in the two- 
body Hilbert space; the biorthogonal partner defined by 
(r\nl) = r -1 (r|nZ). Since the three-body Hilbert space is 
a direct sum of two-body Hilbert spaces, an appropriate 
basis in angular momentum representation is the direct 
product 

\nvl\) a = [\nl) a <gi\v\) a ], (n,u = 0,1,2,...)- (16) 

Here I and A denote the angular momenta associated with 
the coordinates £ and 77, respectively. In this basis the 
completeness relation takes the form 



N 

lim > 



\nvlX) a a (nvlX\ 



lim lft; (17) 

N— >oc 



a sum over the angular momenta being assumed implic- 
itly. Note that in the three-body Hilbert space, three 
equivalent bases belonging to fragmentation a, (3 and 7 
are possible. 

To proceed, it is convenient to approximate Eqs. (|l5| ) 

by 

\4 ±] ) = w<^ )(±) > + GfiU; ] E i&i^). ( 18 ) 

with the short-range potential in the three-body 
Hilbert space taken to have a separable form, viz. 



JV 



\nvlX) a v^ p p{n'v'VX'\ 



(19) 



n.u.n' .v'— 



where v 6 a/3 = a {nvlX\v s a \n'v'l'X') p . In Eq. ((T9|) the ket 
and bra states are defined for different fragmentations, 
depending on the environment of the potential operators 
in the equations. A similar approximation is made on 
the potential U a in Eqs. (jl4|) and (|l^) , with bases of the 
same fragmentation a applied on both sides of the oper- 
ator. Thus, by truncating the short-range operators in 
Eqs. Jl5|), ( |l4| ) and ([l3|), the set of linear integral equa- 
tions reduces to an analogous set of linear algebraic equa- 
tions with the operators replaced by their matrix repre- 
sentations. Calculation of v^p and of U_ a can be made as 



shown previously || . A similar procedure enables calcu- 
lation of the matrix elements G a — a (ni'lX\Ga\n'i/'l'X') a 
with three-body bound states M. The Green's operator 
G Q is a resolvent of the sum of two commuting Hamilto- 



nians, h 



and h T , 



(0 



which act in 



different two-body Hilbert spaces. Thus, using the con- 
volution theorem |fl5|| , which is a direct consequence of 
Dunford's operator calculus |l(| , the three-body Green's 
operator G a equates to a convolution integral of two- 
body Green's operators, i.e. 



G a (z) = (z- he 



2ni 



tt^ f dz' g tac {z - z') g Va (z'), 



(20) 



c 



where g£ a (z) = (z — h^^)^ 1 and g Va (z) = (z — h^ a ) . 
The contour C should be taken counterclockwise around 
the continuous spectrum of h„ a so that g^ a is analytic 
in the domain encircled by C. For bound-state energies 
the spectra of the two Green's operators are well sepa- 
rated and this condition can be fulfilled easily ||. Also, 
below breakup threshold, where the bound-state pole of 
g^ a meets the continuous spectrum of gr) a , contour in- 
tegration can still be performed For other positive 
energy scattering problems, however, the continua over- 
lap so the applied contours are not viable. 

But there exists a contour which is valid for all z of 
physical interest. Besides positive real values of z that 
arise with scattering above the breakup threshold, this 
contour can deal with complex values of z having nega- 
tive imaginary parts that are needed for resonant-state 
calculations. In this approach first one must shift the 
spectrum of g^ a by taking z = E + ie with positive e. 
By doing so, the two spectra are well separated and the 
spectrum of g 71a can be encircled. Next the contour C is 
deformed analytically in such a way that the upper part 
descends to the second Riemann sheet of g Va , while the 
lower part of C can be detoured away from the cut [see 
Fig. [I]] . The contour still encircles the branch cut singu- 
larity of g Va , but in the e — > limit it now avoids the 
singularities of g^ a . Moreover, by continuing to negative 
values of s, the branch cut and pole singularities of g^ a 
move onto the second Riemann sheet of g Va and, at the 
same time, the branch cut of g^^ moves onto the second 
Riemann sheet of g^ a . Thus, the mathematical condi- 
tions for the contour integral representation of G a (z) in 
Eq. (|2^) can be fulfilled for all energies. In this respect 
there is only a gradual difference between the bound- and 
resonant-state calculations and scattering ones below and 
above the breakup threshold. 

The matrix elements G ' can be cast in the form 



G a (z) = — f dz' a (nl\g (a (z ~ z')\n'l') a 



c 



x a{v>\9T, a {z')\v'\') a . 



(21) 
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The matrix elements of the two-body Green's functions 
in the integrand are known for all complex energies from 
two-body state properties jl7],[l8) . In this formalism the 
Faddeev components appear as linear combinations of 
functions 

(^aVa\tpa) ~ ((.aVa \ G a (z) \nvl X) a = 

-^j> dz' (£ a \g ia (z - z')\nl) a {■q a \g Va {z')\u\) a , (22) 
which are convolution integrals of Coulomb- like functions 




In this paper I have shown that the set of Faddeev 
and Lippmann-Schwinger equations for the three-body 
Coulomb problem derived from the "three-potential" pic- 
ture H possess compact kernels for all energies. I have 
found an analytic representation of the channel-distorted 
Coulomb Green's operators in terms of a convolution in- 
tegral of two-body Green's operators. The method facil- 
itates the solution of integral equations in which the ex- 
act bound- scattering- or resonant-state Coulomb asymp- 
totics are automatically incorporated. Solution of these 
equations in a Coulomb-Sturmian space representation 
is most practical, as then there is an analytic representa- 
tion of the terms in the contour integrals. 
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FIG. 1. Analytic structure of g^ a (z — z') g Va (z') as a func- 
tion of z with z = E + ie, E > 0, e > 0. The contour C 
encircles the continuous spectrum of h Va ■ A part of it, which 
goes on the second Riemann-sheet of h Vct , is drawn by broken 
line. 
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